Density Matrix Functional Theory for the Lipkin model 
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A Density Matrix Functional theory is constructed semi-empirically for the two-level Lipkin model. 
This theory, based on natural orbitals and occupation numbers, is shown to provide a good descrip- 
tion for the ground state energy of the system as the two-body interaction and particle number vary. 
The application of Density Matrix Functional theory to the Lipkin model illustrates that it could 
be a valuable tool for systems presenting a shape phase-transition such as nuclei. The improvement 
of one-body observables description as well as the interest for Energy Density Functional theory are 
discussed. 
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I. INTRODUCTION 

Recently, large efforts are devoted to the construc- 
tion of an Energy Density Functional (EDF) able to de- 
scribed at best properties of nuclei over the whole nuclear 
chart [E [3. The standard strategy to design an EDF 
for nuclei is to start with a single-reference EDF (SR- 
EDF) where an effective interaction (Skyrme or Gogny 
type) and a trial state (Slater Determinant or more gen- 
erally quasi-particle vacuum) are chosen. This tech- 
nique is able to describe short-range correlations like 
pairing and provides already a rather good description 
of observables such as masses, under the condition that 
some symmetries of the original Hamiltonian are bro- 
ken. The SR-EDF is then extended to restore broken 
symmetries and/or incorporate long-range correlations 
through configuration mixing, leading to the so-called 
Multi-Reference EDF (AIR-EDF) 0] . Recent applications 
of this technique have revealed important conceptual and 
practical difficulties Q related to the absence of a con- 
structive framework for multi-reference calculations. A 
solution to this problem, has been recently proposed @ 
and successfully tested in nuclei However, this cure 
does not apply to most of the functionals currently used 
01 , i.e. those with fractional powers of density. 

This motivates the search of new techniques to extend 
actual SR-EDF. The Density Matrix Functional Theory 
(DMFT) Q appear as an alternative to configuration 
mixing [9| . Although this theory was proposed more than 
30 years ago Q, explicit forms of functionals and appli- 
cations have only been explored rather recently. There 
are nowadays an increasing interest in proposing accurate 
DMFT In this work, DMFT is applied to the two- 
level Lipkin model In this model, the Hartree-Fock 
(HF) theory fails to reproduce the ground state energy 
|l2l | while configuration mixing like Generator Coordi- 
nate Method (GCM) provides a suitable tool 0]. There- 
fore, the two-level Lipkin model is perfectly suited both 



to illustrate that DMFT could be a valuable tool and 
to provide an example of functional for system with a 
"shape" phase-transition. 

In this following, some aspects of DMFT are first re- 
called. Then a semi-empirical functional is constructed 
for the Lipkin model and applied for various particle 
number and two-body interaction strengths. It is shown 
to improve significantly the HF theory. Finally, the in- 
terest of constructing more general functional of natural 
orbitals and occupation numbers in the EDF context is 
outlined. 



II. DISCUSSION ON DMFT 

The concept of Density Matrix Functional Th eory is a 
generalization of the Hohenberg-Kohn theorem [IJl due 
to Gilbert Q. It relies on a theorem showing that the 
ground state energy could be written as a functional of 
the one-body density matrix (OBDM) 7(r, r') (instead of 
the local one-particle density p{r) = 7(r, r) in the stan- 
dard Hohenbcrg-Kohn theorem). Then, similarly to the 
Kohn-Sham orbitals [l^ . the eigenvalues n.i and eigen- 
vectors ifi, called hereafter resp. occupation numbers 
and natural orbitals, of the OBDM are often used in- 
stead of 7(r, r') with the relation 7 = J2i \'fi)n'i{'Pi\- The 
variation of the functional 



- ^i{Tr{p) - N} 
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with respect to one-particle state components (p* (r) and 
occupation numbers (with the additional constraint < 
ni < 1) is then performed to obtained the optimal (pi, 
rii and associated ground state energy. The set of La- 
grange multipliers p, and {A.^} are introduced to insure 
particle number conservation and orthogonality of the 
single-particle states. In Eq. ([1]), £[{</Ji}, {n^}] is nothing 
but the functional itself which has to be found. In elec- 
tronic system, the functional is generally separated into 
the Hartree, denoted by Sh, (eventually Hartree-Fock, 



2 



£hf) part and the exchange-correlation part, denoted 
here by £xc (eventually correlation only Ec)- 

While the DMFT has been studied theoretically for a 
rather long time [1, [l^, [l^, [3 ; only recently explicit 
functionals of the OBDM or directly to natural orbitals 
have been proposed and applied to realistic situations [13, 
M: HO, HH, M: In, 113, m, m, M, Hi, Hi, H [MI . There are 
nowadays extensive works to test functionals especially in 
infinite systems, the so-called Homogeneous Electronic 
Gas (HEG) 



III. APPLICATION OF THE DMFT TO THE 
LIPKIN MODEL 

The "Lipkin Model" [ll| is an exactly solvable model 
that has often been used as a benchmark for approxi- 
mations for the nuclear many-body problem Q. In this 
model, one considers N particles distributed in two N- 
fold degenerated shells separated by an energy e. The 
associated Hamiltonian is given by: 



H = e.h-^{J+J+ + J-J^), 



(2) 



where V denotes the interaction strength while Jq, J± 
are the quasi-spin operators defined as 



1 ^ 



0=1 



N 



p and cL p are creation operators associated with the 
upper and lower level. The exact solution of this model, 
is easily obtained by noting that (but not Jq) commute 
with H. It is then convenient to introduce the basis of 
eigenstates of and Jq and diagonalize the Hamiltonian 
in this particular space (for more detail see for instance 



automatically insures the orthogonality of the new states. 
Due to simple structure of the Lipkin model, the variation 
with respect to the SD state is identical to the variation 
of the a and (p, i.e. |$) = \^{a,ip)) and the HE energy 
becomes a functional of these parameters: 



(5) 



Anticipating for the forthcoming discussion, we first write 
£mf as a functional of the OBDM 7: 



SmfIj] = eTr(Jo7) 
V{N - 1) 
27V 



{(Tr[7J+])2 + (Tr[7J-])2}.(6) 



In the Hartrce-Eock limit, the OBDM contains all the 
information on the many-body state and simply reads 

7 = Y.p=i \fo,p){fo,p\- Reporting the expression of |cpo,p) 
in terms of a and tp, we recover the standard HE expres- 
sion: 



£mf [a, (fi] 



zN 



|cos(2a) 4- I sm^{2a) cos(2(^)|(.7) 



where x = V{N—l)/e. Minimizing with respect to (a, ip) 
leads to HE energy, denoted by £^p (both with = 0): 



HF 



sN 
sN 



for X < 1 (at a = 0), 

,2^ 



£hf = (1 + X') for X > 1 (at xcos(2a) = 1). 

The HE solution for the Lipkin model has been exten- 
sively discussed in the literature While it provides 
a rather good estimate of the exact energy in the weak 
coupling or in the large N limit, it generally differs rather 
significantly from it in particular for x ^ 1- This dis- 
crepancy essentially reflects the failure of the HE method 
to account for configuration mixing in a single-reference 
framework for systems with a shape phase-transition . 



B. Expression of the Hamiltonian for a General 
correlated state 



A. Hartree-Fock approximation 

In the Hartree-Eock (or Mean-Eield) theory, the many- 
body wave function is replaced by a Slater Determi- 
nant (SD) given by |$) = n^^ajp]— ). Here, a new 
single-particle basis, denoted by {po,p, Vi,p} associated 
to the set of creation/annihilation operators {«op'^ip} 
has been introduced through the relation 




where the choi 



Due to the two-body nature of the Hamiltonian (Eq. 
(21)), the most natural way to extend the Mean-Eield 
framework to correlated system is to introduce the two- 
body density matrix, denoted by F12 and the associated 
correlation matrix ai2 (see for instance [35|). Using the 
OBDM of the correlated system 7, (T12 is defined through 
the relation: 

= 7172(1 - fl2) + 0-12, (8) 

where P12 denotes the anti-symmctrization operator 
while the label "i" in ji refers to the particle on which 
the density is applied, i.e. {ij\jij2\kl) = (*|7|fc)(j|7lO 
[36t . The expectation value of the energy then splits into 
a mean-field part and a correlated part as : 



/ = cos(a), g = sin(Q;)e*'^, 



(4) 



(9) 
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where ^mfM is given by expression ([6]) except that 7 
now refers to the OBDM of the correlated system. £c 
reads: 

V 

^c[^12] = -yTr{[J+ J+ + J_]i2Cri2} 



where the notation [.] 12 emphasizes that the trace is per- 
formed using the two-body matrix elements of the oper- 
ators. 

Eq. ([9]) is valid for any correlated system including the 
exact ground state. The minimization of it with respect 
to any kind of OBDM and correlations will therefore lead 
to the ground state energy. Such a direct minimization is 
complex due to the number of degrees of freedom involved 
in the variation [3^. 

DMFT provides a practical solution to this problem. 
Indeed, according to this theory Q, the total energy 
could be written as a functional of 7. Since £mf[i\ is 
already written as a functional of the OBDM, we are left 
to find a functional for £c- In many cases, both £mf 
and £c are directly written as a functional of the natural 
orbitals (fi and occupation numbers Ui, i.e. 

£[{(Pi,ni}] = £MF[Wi,ni}]+ £c[{'Pi,n.i\]. (11) 

DMFT has some additional advantages compared to 
standard DFT. Besides the fact that the exchange con- 
tribution could be expressed exactly with the OBDM, 
at the minimum the optimal OBDM identifies with the 
ground state OBDM. Therefore, expectation values of 
any one-body observable could be computed and should 
correspond to the ground state expectation value. Ac- 
cordingly, if the mean-field prescription is used in £mf 
for a given H, the value of £mf at the minimum will 



truly correspond to the true mean-field contribution to 
the total energy. Consequently, £c will truly correspond 
to the contribution of correlation (this issue will further 
be discussed in section UlI Ep . 

Recently, significant efforts have been made in the 
practical construction of such density matrix functionals 

m E s mi, m, m, h m, m, ^ m, m, im. in 

most cases, guided by general considerations [18[ and/or 
constraints on the two-body density [IBl , cri2 is first writ- 
ten as a functional of the natural orbitals and occupation 
numbers. This corresponds generally to a starting point 
for more elaborated functionals. Functionals are further 
enriched by incorporating additional terms either to cor- 
rect for the self-interaction problem or eventually to bet- 
ter reproduce some specific limits in the infinite HEG or 
configuration interaction calculation in molecules (for a 
recent review see [3?}). 

C. Construction of a density matrix functional for 
Correlated state 

Due to the specific form of the Lipkin Hamiltonian 
given by Eq. ([2]), 7 simply writes in the natural basis as: 



N 

7 = ^{\'PQ,p)'n'Q{vo.p\ + |v5i,p)ni((/?i,p||, (12) 
p=i 

with ni = (1 — no). The single-particle states \(pi,p) (with 
i — 0, 1) now stand for the natural orbitals while the 
correspond to occupation numbers. Similarly to the HF 
theory, creation operators, denoted by a|p, associated to 

natural orbitals are expressed from the c^j_ ^ using Eq. 
([3]). The mean-field contribution is easily deduced from 
the HF case, using Eq. ^ and the expression of 7 given 
above, leads to: 



£MF{{'-Pi,p,ni}) 



-|7V| cos(2a)(2no - 1) + | sin2(2a) cos{2tp){2na - 1)2|. 



(13) 



Expressing £c is less straightforward. A possible strat- 
egy to construct functionals is to identify specific limits 
at which explicit forms are known. For instance the two 
electron case which has been studied in has largely 
influence presently used DMFT for molecules [IJ- Fol- 
lowing this idea, the N = 2 case is first considered. 

1. The N = 2 case 

To study the N = 2 case, a basis of Slater Determinant 
is constructed using the natural orbitals, i.e. = 



o-i pttj with i and j cither or 1. The ground state 
wave-function ^E* then reads: 

|*)=5]C,,|<i>,;,). (14) 

From the expression of the OBDM and using the fact, 
that the single-particle states are natural orbitals we ob- 
tain the simple relation |CyP = SijUi. From which we 
deduce Cij = Sije^^" ^/ui. As illustrated below, the sim- 
plest choice e**" = 1 is convenient and leads to 

I*) = V^|$oo) + (15) 



4 

which is nothing but the exact ground state wave- leads to a total ground state energy £ given by: 
function written as a functional of occupation numbers 
and natural states. Reporting this expression in ($|lf|$), 



£ (a, if, no) = -£cos(2a)(2no - 1) - vj- sin2(2a) cos(2(/3) + 2 (sin''(a) cos(4(p) + cos'^(a)) V?io(l - no)j{16) 
Using the expression of the mean- field contribution (Eq. (fT3|) ). we deduce 

£q {a, if, no) — — 2f| sin^(2a) cos(2i^)no(l — no) + (sin''(Q!) cos(4iy9) -f- cos''(a)) y/ no{l — no)|. (17) 

I 




FIG. 1: Comparison between the exact energy (solid line), 
Hartree-Fock (dashed line) and the energy obtained from the 
minimization of Eq. (|16|l (filled circles) as a function of x for 
the N = 2 case. 

Since the above functional is exact, the ground state en- 
ergy should be recovered by minimization with respect to 
no, a and ip. As in the HF case, cp = could always be 
taken. The variation of no should be made under the con- 
straint no € [0, 1]. A similar technique as in the HF-I-BCS 
case could be employed Writing no — cos^{6) with 
9 G [0, 7r/2]. gives at the minimum 

^ ^ / sin^(a)cos(4<,) +cos^(a)\ 
\ C0s(2q;) / 

Then, variation of the a is made to obtain the minimum 
energy. The result (filled circles) is displayed in Fig. [T] 
and compared to the exact ground state energy given by 
E ~ ~ \/l + [m (solid line) and the Hartrce-Fock en- 
ergy (dashed line). Not surprisingly, while the HF curve 
deviates significantly from the solid line, the DMFT could 
not be distinguished from the exact case. 

2. DMFT for N > 3 and large N scaling 

The Lipkin model with TV = 2 is an interesting peda- 
gogical example of DMFT where the exact energy func- 



tional in terms of natural orbitals and occupation num- 
bers is known. This limit is used here as a guide to 
provide a DMFT for > 3. The most simple extension 
of the functional derived for N = 2, consists in assuming 
that the interaction of the A^ particle could be written 
as a sum of the interaction of the A^(A^ — l)/2 pairs of 
particles, each pair contributing to the total energy as in 
the N = 2 case. This prescription naturally leads to the 
mean-field Hamiltonian given (fTH]) and amount to take 
for all N 

, A^(A^-l) «=2, 
£c(a,^p,no) = £c (a,(p,no). (19) 

The minimal total energy obtained by varying both oc- 
cupation numbers and a (using ip = 0) using this pre- 
scription is displayed in Fig [2] (solid circles) as a func- 
tion of X foi' various particle numbers. In each case, the 
exact solution (solid line) and the Hartree-Fock prescrip- 
tion (dashed line) are shown. The simple scheme using 
(|19p clearly provides a very poor approximation for the 
ground state energy and always leads to an energy much 
below the exact one. In addition, the discrepancy in- 
creases as A^ increases. This failure points out the com- 
plex many-body correlations present in the Lipkin model 
coming from the mixing of 1 particle-1 hole (Ip-lh), 2p- 
2h, np-nh excitations in the component in the ground 
state. This leads to a much more complex situation than 
the N ^ 2 case (Eq. (jTSj) ). In particular, using the con- 
tracted Schrodinger equation (CSE), we do expect that 
the two-body correlation depends on higher order corre- 
lation matrices [22j . The correlation energy given by (jl9p 
clearly neglects these higher-order effects. 

From Fig. [21 we see that the correlation energy is 
largely overestimated. In the following, we use a slightly 
different prescription for the correlation energy given by: 

N>3 A^(A^ — 1) JV = 2 

£c ia,(p,no) = ri{N) (a,(p,no), (20) 

where r](N) is a reduction factor introduced to mimic the 
effect of higher order correlations. The optimal value of 
r] this is determined empirically from the following pro- 
cedure. For a given value of A^ and 77, the mimimum en- 
ergy of the corresponding DMFT, denoted by £,^j„(ry, x) 
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FIG. 3: Values of the optimal quenching factor rj{N) as a 
function of the particle number N. The solid line represents 
the function ri{N) = c Af"^/^ with c = 1.5. 



creases. Since this correlation energy is proportional to 
(J^), we observe that the N~^^^ deduced empirically is 
nothing but the la,rge N scaling behavior obtained ana- 
lytically in ref. lia], i.e. ( j2)/iV2 oc 7V-2/3. 



FIG. 2: Exact ground state energy (solid lines) displayed as 
a function of x for A'^ = 5 to 20 resp. from top to bottom. 
In each case, the corresponding HF (dashed line) and DMFT 
(filled circles) are shown. The latter are obtained by mini- 
mization of the functional using the mean-field and correla- 
tion energy resp. given by Eq. (|13p and Eq. (|19p . 



is computed for x between and Xmax — 3. Then, the 
quantity D^{t]), given by 

D'^iv) - {£^six)-£lniv,x)Vdx, (21) 

Jo 

where £Qg{x) denotes the exact ground state energy 
for a given N and Xi is computed. Obviously, {rf) 
gives a measure of the deviation between the DFMT 
minimum energy and the exact energy over the inter- 
val X S [0, Xmaa:]. Thc Optimal value of ri{N) is defined 
as the minimum of (jj) as 77 varies between (the HF 
limit) and 1 (the prescription of Eq. (jTO]) ). Values of 
optimal reduction factors are reported by filled circles in 
Fig. [3] as a function of A''. As guessed from the increasing 
discrepancy observed with increasing TV observed in Fig. 

the higher is A'^ the smaller should -q be taken. Thc 
variation of "qiN) turns out to simply behave as N^^^^. 
The solid line in Fig ^ represents the function 

ri{N) = c iV"2/^ (22) 

with c = 1.5 deduced by fitting optimal values of ?/. In 
the following, we show that the semi-empirical density 
matrix functional theory constructed from the mean-field 
and correlation energy resp. given by Eq. (|13[) and Eq. 
(PO)) significantly improves the HF theory. 

Combining expression with Eq. (PO)) shows that 
the correlation energy scales as £c oc N^^"^ as N in- 



D. Results of the semi-empirical DMFT 

Thc minimal energy deduced from the semi-empirical 
density matrix functional proposed above (filled circles) 
is systematically compared to the exact ground state en- 
ergy (solid line) and HF energy (dashed line) for differ- 
ent particle number and two-body interaction strength 
in Fig. |4]and[5l In all cases, the DMFT significantly im- 
proves the HF result and turns out to be very close to the 
exact one. As illustrated in Fig. [4l the HF energy gen- 
erally deviates rather significantly from the exact energy 
around x = 1 ^^id does not provides the correct asymp- 
totic behavior as x increases. This deviation, which dis- 
appears as TV —> 00, is due to the failure of HF theory 
in thc presence of a "shape" phase transition [l3| be- 
tween the spherical solution (x < 1) and the "deformed" 
solution (x > 1). The standard technique to properly 
account for this effect is to mix different slater determi- 
nants like in the GCM theory. We see in figure [5] and [5] 
that, except a small deviation around x = 1 which seems 
to slightly increase as TV increases, both the asymptotic 
behavior at large x and TV and the energy around x = 1 
arc rather well reproduced. It is worth mentioning that 
the DMFT is much less demanding in terms of computa- 
tional power than thc GCM and thus provides a rather 
interesting alternative to the latter theory. 



E. Discussion on Density Matrix Functional 
Theory with broken symmetry 

Similarly to the Hartrce-Fock case, the one-body den- 
sity solution of the functional developed here can eventu- 
ally violate some symmetries of the "true" ground state 
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FIG. 4: Exact ground state energy (solid lines) displayed as 
a function of x for A'^ = 5 to 20 resp. from top to bottom. 
In each case, the corresponding HF (dashed line) and DMFT 
(filled circle) minimum energy are shown. The DMFT calcula- 
tion is performed using the mean-field and correlation energy 
resp. given by Eq. ((T3} and Eq. ^ with t){N) = 1.5 A^"^^^. 



FIG. 5; Exact ground state energy per particles (solid lines) 
displayed as a function of A'^ for x ~ '2 and 3 from top to 
bottom. In each case, the corresponding HF (dashed line) and 
DMFT (filled circle) minimum energy are shown. The DMFT 
calculation is performed using the mean-field and correlation 
energy resp. given by Eq. (|13|l and Eq. (|2Up with 77(7V) = 
1.5 N-'^/^. 



density. The invariance of the Hamiltonian ([2]) with re- 
spect to parity imposes (J+) = {J-) = 0. This imphes 
that at the minimum of the functional a = 0. This is in- 
deed the case for the exact functional given by Eq. (fTBl) 
for N — 2. However, solutions with a ^ are found 
for larger particle number. This is illustrated in figure [5] 

where the value of A_| = (( J+) -I- {J-))/2N is displayed 

as a function of x- The value of x, for which a becomes 
non-zero, denoted by Xc is infinite for N = 2, around 1.6 
for N = 5 and tends to the HF case (xc = 1) as goes to 
infinity. Note that, in this limit, the HF functional alone 
already provides a very good functional for the Lipkin 
model. 

According to DMFT, the OBDM 7 obtained by min- 
imizing the exact functional should match the exact 
OBDM at the minimum. Therefore, strictly speaking, 
the extracted one-body density and associated natural 
states and occupation numbers could not be the exact 
OBDM when some of the symmetries of the system are 
broken. This has to be kept in mind when discussing 
expectation values of one-body observables. As an illus- 
tration, the expectation value of the single-particle part 
of the Hamiltonian £j„ = eTr{'yJo) obtained using the 
OBDM minimizing the DMFT is displayed in Fig. [7] as a 
function of x for different particle number (filled circles). 
The exact (solid line) and HF (open circles) prescription 
are also shown. The value of the total energy minus Sj^ 



is also shown in each panel for the DMFT (filled square) , 
exact (dashed line) and HF (open square). 

The density obtained in DMFT using the semi- 
empirical functional described in section fill CI always im- 
proves the estimate of £j„ for small x- It generally gives 
an almost perfect result when a = at the minimum, i.e. 
when the parity symmetry is respected by the OBDM. 
The behavior of Ej^ is also rather satisfactory at large N 
and X- In other cases, small particle number and large x 
or large particle number and intermediate X (1 ^ X ^ 2). 
Fig. [7] illustrates that, when the functional respects all 
symmetries of the original Hamiltonian, expectation val- 
ues of one-body observables perfectly match the exact 
results. This is the case of the semi-empirical functional 
proposed here for all particle number and x ^ 1- Of 
course it would be desirable to provide functionals that 
respects the symmetries at the first place. However, as it 
is well known already at the HF level, the introduction 
of theories where symmetries are explicitly broken is a 
way to grasp some of the correlations which would have 
been very hard to incorporate without breaking the sym- 
metries. The success of the DMFT functional at large N 
and X can be attributed to the explicit parity symmetry 
breaking as in the HF case. 
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FIG. 6: Expectation value of A-| as a function of x for N = 5 

(top), iV = 10 (middle) and iV = 20 (bottom). The DMFT 
result (filled circles) is systematically compared with the HP 
(open circles) value. 



IV. SUMMARY AND DISCUSSION ON EDF 

In this work, the Density Matrix Functional Theory 
is applied to the two-level Lipkin model. Guided by the 
N ^ 2 case, a semi-empirical functional of the natural 
wave-functions and occupation numbers is constructed. 
The minimization of the DMFT is shown to give a much 
better agreement of the exact ground state energy than 
the HF scheme over a wide range of particle number and 
two-body interaction strength. The success of DMFT 
in the Lipkin model, shows that this theory could be a 
valuable tool in many-body system in the presence of 
"shape" phase transition. Such transition often occurs 
for instance in nuclear physics [l|, Q and is generally 
treated by first introducing the Energy Density Func- 
tional of a single-reference vacuum (Slater Determinant 
or quasi-particle state) and then use GCM theory Q. 

DMFT could be a powerful tool to improve actual SR- 
EDF functionals, by writing directly the EDF in terms 
of natural orbitals and occupation numbers. The pos- 
sibility to introduce occupation numbers has been pro- 
moted in Ref. [i^l for the pairing Hamiltonian and in 
for the three-level Lipkin model using a slightly dif- 
ferent technique. However, the strategy based on DMFT 
seems quite natural to extend actual EDF. Indeed, fol- 
lowing the strategy used here for the Lipkin model, we 
can write the nuclear EDF as 

£EDF[Wi,n^}] ^ £MF[{Vi,ni}] + £c[{Vi,nt}]- (23) 



FIG. 7: Expectation value of the one-body part of the Hamil- 
tonian, denoted by £jg obtained at the minimum of the 
DMFT as a function of x for = 5 (top), iV = 10 (mid- 
dle) and iV = 20 (bottom). The DMFT result (filled circles) 
is systematically compared with the exact (solid line) and HF 
(open circles) value. The value of the total energy minus £jg is 
also shown in each panel for the DMFT (filled square), exact 
(dashed line) and HF (open square). 

The most natural choice for the mean-field part is to use 
the actual Skyrme functional which has been optimized 
for decades. Note that, the nuclear problem differs from 
the electronic case and/or the Lipkin model presented 
here due to the fact that coefficients of the functional are 
not directly linked to the bare interaction but adjusted 
on experimental data. Therefore, the mean-field contri- 
bution already contains a large fraction of the correla- 
tion. Nevertheless, the above decomposition (Eq. ([^5]) ) 
is already used in the nuclear context in SR-EDF cal- 
culations when a quasi-particle vacuum is retained for 
the trial state. Then, the correlation energy identifies 
with the pairing energy which could, in the canonical 
basis, be written as a functional of occupation numbers 
and natural orbitals. Indeed, using the notation (i,i) 
for canonicals pairs of single-particle states, the pairing 
energy writes 

(24) 

where v'^'^ is the effective interaction in the pairing 
channels and where we have replaced components of 
the anomalous density k in the natural basis by = 
•ni{\ — Ui). This functional is adapted for pairing like 
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correlations. However, the use of a reference state writ- 
ten as a product of quasi-particle states clearly restricts 
the type of density matrix functional that can be guessed. 
This class appears to be not general enough to account 
for the diversity of phenomena occurring in nuclei. The 
possibility to use different functional than the BCS like 
ones has already been discussed in the early time of the 
Skyrmc EDF history [i^ . The functional developed here 
for the Lipkin model as well as functionals recently pro- 
posed in electronic systems [131, clearly point out the 
possibility to use alternative functionals which could be 
of interest for nuclear systems. 

A crucial aspect of the present work is the introduc- 
tion of DMFT functionals that explicitly breaks some of 
the symmetries of the original Hamiltonian to incorpo- 



rate complex correlations. It should be kept in mind that 
broken symmetries imply that symmetries should a pri- 
ori be restored. The problem of restoration of broken 
symmetries in functional theories is an important aspect 
which deserves specific studies in the near future [5|, l6|, |3] ■ 
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